arXiv:1509.03345v2 [hep-ph] 14 Sep 2015 


The magnetized effective QCD phase diagram 


Alejandro Ayala^’^, C. A. Dominguez^, L. A. Hernandez^, M. Loewe^’^’^, R. Zamora^ 
^Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 

Apartado Postal 70-543, Mexico Distrito Federal 04510, Mexico. 

^ Centre for Theoretical and Mathematical Physics, and Department of Physics, 
University of Cape Town, Rondebosch 7700, South Africa 
^Institute de Fisica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile. 
'^Centro Cientifico-Tecnologico de Valparaso, Casilla 110-V, Valparaiso, Chile. 

The QCD phase diagram in the temperature versus quark chemical potential plane is studied in 
the presence of a magnetic field, using the linear sigma model coupled to quarks. It is shown that the 
decrease of the couplings with increasing field strength obtained in this model leads to the critical 
temperature for the phase transition to decrease with increasing field intensity (inverse magnetic 
catalysis). This happens provided that plasma screening is properly accounted for. It is also found 
that with increasing field strength the location of the critical end point (CEP) in the phase diagram 
moves toward lower values of the critical quark chemical potential and larger values of the critical 
temperature. In addition, the CEP approaches the temperature axis for large values of the magnetic 
field. We argue that a similar behavior is to be expected in QCD, since the physical impact of the 
magnetic field, regardless of strength, is to produce a spatial dimension reduction, whereby virtual 
quark-antiquark pairs are closer on average and thus, the strength of their interaction decreases due 
to asymptotic freedom. 
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I. INTRODUCTION 

The behavior of strongly interacting matter in presence 
of magnetic fields has become a subject of increasing in¬ 
terest over the past few years. This interest has been 
sparked by recent lattice QCD (LQCD) results showing 
that the transition temperature with 2-1-1 quark fla¬ 
vors decreases with increasing magnetic field and that 
the strength of the condensate decreases for tempera¬ 
tures close and above the phase transition EHi. This be¬ 
havior has been dubbed inverse magnetic catalysis, and 
it has been the center of attention in a large number 
of mo del- dependent analyses [Ml. For recent reviews 
see Refs. [H HI. In general terms, it seems that in¬ 
verse magnetic catalysis is not obtained in mean field 
approaches describing the thermal environment [IMl, 
nor when calculations beyond mean field do not include 
magnetic effects on the coupling constants HI. 

The novel feature implemented in effective models, 
able to account self-consistently for inverse magnetic 
catalysis, is the decrease of the coupling constants with 
increasing field strength obtained from the model it¬ 
self without resorting to ad hoc parametrizations. 

This has been achieved within the Abelian Higgs model 
and the linear sigma model coupled to quarks (LSMq). 
This behavior is made possible by accounting for the 
screening properties of the plasma, which have been re¬ 
cently formulated consistently for theories with sponta¬ 
neous symmetry breaking |25| . This results in a formal¬ 
ism beyond the mean field approximation [2^ . Screen¬ 
ing is also important to obtain a decrease of the coupling 
constant with the magnetic field strength in QCD in the 
Hard Thermal Loop approximation |28j . It has also been 


shown, by means of a QCD sum rules analysis at zero 
temperature, that both, the threshold energy for the on¬ 
set of the continuum in the quark vector current spec¬ 
tral density (a phenomenological parameter that signals 
the onset of deconfinement) and the gluon condensate 
increase with increasing magnetic field, as expected . 

Recently, the LSMq has been also used to explore the 
phase diagram without magnetic fields [s^ . It was found 
that there are values for the model couplings that allow 
locating a critical end point (CEP) in the region where 
lattice inspired calculations find it M- Since the LSMq 
does not exhibit confinement, this behavior is attributed 
to the proper treatment of plasma screening, instead of 
to the existence of a given confinement length scale [s^. 

A pertinent question is whether the above description 
in the presence of a magnetic field can be used to study 
how such CEP changes with the field intensity. Recent 
LQCD calculations show that for very strong mag¬ 
netic fields, inverse magnetic catalysis prevails and the 
phase transition becomes first order at asymptotically 
large values of the magnetic field for vanishing quark 
chemical potential yi. A similar behavior is obtained in 
the Nambu Jona-Lasinio model if one includes a magnetic 
field dependence of the critical temperature in agreement 
with LQCD [s^. 

In this paper we use the LSMq to explore the con¬ 
sequences of a proper handling of the plasma screen¬ 
ing properties in the description of the magnetized ef¬ 
fective QCD phase diagram. We show that when includ¬ 
ing self-consistently magnetic field effects in the calcu¬ 
lation of both the effective potential as well as on the 
thermo-magnetic dependence of the coupling constants, 
the CEP’s location moves toward smaller values of the 
critical quark chemical potential, and larger values of the 
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critical temperature. In addition, above a certain value 
of the field strength the CEP moves to towards the T- 
axis. We argue that this behavior can be understood on 
general grounds, as the magnetic field produces a dimen¬ 
sion reduction, whereby virtual charged particles from 
the vacuum are effectively constrained to occupy Lan¬ 
dau levels, thus restricting their motion to a plane. This 
makes these particles to lay closer to each other on aver¬ 
age, thus reducing the interaction strength for strongly 
coupled theories. This situation takes place regardless of 
how weak the external field is. 

The paper is organized as follows: In Sec. HI] we recall 
the basic features of the linear sigma model with quarks, 
and in the presence of a magnetic field. Next, we write 
the effective potential at finite temperature and quark 
chemical potential, including all the degrees of freedom 
of the model. In Sec. uni we find the thermo-magnetic 
corrections to the boson self-coupling and to the cou¬ 
pling between fermions and bosons. In Sec. IIVI we find 
the effects of the thermo-magnetic dependence of the cou¬ 
plings on the critical temperature for chiral symmetry 
restoration transition. We find that Tc is a decreasing 
function of the magnetic field for any value of /i. We also 
study the magnetized effective QCD phase diagram and 
in particular find how the CEP’s location changes with 
the field strength. We finally summarize and conclude in 
Sec.lVj In the appendix we describe the computation of 
the thermo-magnetic corrections to the couplings, and to 
the temperature and density dependence of the fermion 
thermal mass. 


II. EFFECTIVE POTENTIAL 

The Lagrangian of the sigma model, including quark 
degrees of freedom, is given by 




-I- — g'tp{a + i^ 5 f ■ Tr)tp, 


( 1 ) 


where ip is an SU(2) isospin doublet, tt = (tti, 7r2 ,tts) is 
an isospin triplet and a is an isospin singlet, with 




( 2 ) 


is the covariant derivative. is the vector potential cor¬ 
responding to an external magnetic field directed along 
the z axis. In the symmetric gauge it is given by 


= y(0,-y,a:,0). 


(3) 


where q is the particle’s electric charge. satisfies the 
gauge condition = 0. The gauge field couples onlj 
to the charged pion combinations, namely 


7r± = -^ (tti =F W2). 


(4; 


The neutral pion is taken as the third component of 
the pion isovector, 7r° = The gauge field is taken 
as classical and thus we do not consider loops involving 
the propagator of the gauge field in internal lines. The 
squared mass parameter and the self-coupling A and 
g are taken to be positive. 

To allow for spontaneous symmetry breaking, we let 
the a field to develop a vacuum expectation value v 

a ^ a + V, (5) 


This vacuum expectation value can later be identified as 
the order parameter of the theory. After this shift, the 
Lagrangian can be rewritten as 

^ ^ (3Au^ - a^) 

- + iqA^Y^] - i - a^) -h 

- + iipj'^Df_,ip - gvipip + (6) 


where and £j are given by 
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The terms given in Eq. ([7]) describe the interactions 
among the fields a, 7? and ip, after symmetry breaking. 
Erom Eq. (jG]) one can see that the a, the three pions and 
the quarks have masses given, respectively, by 


= 3Au^ — a^, 

= \v —a , 

TUf = gv. (8) 


Using Schwinger’s proper-time method, the expression 
for the one-loop effective potential for one boson field 
with squared mass and absolute value of its charge 
Qb at finite temperature T in the presence of a constant 
magnetic field can be written as 




d°k r°° ds 
(27r)3 _/o cosA{qbBs) 




(9) 


where ujn = 2n7iT are boson Matsubara frequencies. 
Similarly, the expression for the one-loop effective po¬ 
tential for one fermion field with mass m/ and absolute 
value of its charge qj at finite temperature T and chem¬ 
ical potential g, in the presence of a constant magnetic 
field can be written as 





dTV?f 


d°k r°° ds 
(27r)3 Jq cosh(qfBs) 


X 




tanh(gj? Bs) 


+m^j+rqfB] 


(10) 
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where = (2n + 1)ttT are fermion Matsubara frequen¬ 
cies. The sum over the index r corresponds to the two 
possible spin orientations along the magnetic field direc¬ 
tion. 

Including the u-independent terms, choosing the renor¬ 
malization scale as and after mass and charge 

renormalization, the thermo-magnetic effective potential 
in the small to intermediate field regime, in a high tem¬ 
perature expansion can be written as 
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where (x) is the digamma function, Lz„ is the polylog¬ 
arithm function of order n, q is the absolute value of the 


charged pion charge (g = 1), qu = 2/3, and q^ = 1/3 are 
the absolute values of the up- and down-quarks, respec¬ 
tively, 7 b is the Euler’s gamma, Nf = 2 is the number of 
light-quark flavors, and iVc = 3 is the number of colors. 
In order to obtain the leading magnetic field contribution 
of Eq. (ITT]) we use the Euler-MacLaurin expansion [2^ 
in Eqs. dH) and dlo}. 

Though we take the quark masses as equal, the nota¬ 
tion emphasizes that the effective potential is evaluated 
after accounting for the different quark charges. We have 
introduced the leading temperature plasma screening ef¬ 
fects for the boson masses squared, encoded in the boson 
self-energy 11. The leading contribution for the boson 
self-energy in a high temperature expansion, and at fi¬ 
nite fi is given by 

n = aE _ 7VyW<7'E[L*2(-e^/^) + Lz2(-e-'^/^)]. 

Z TT 

( 12 ) 

For the Hurwitz zeta function ({—l/2,z) in Eq. dTlT) to 
be real, one needs 


— -b n > qB, (13) 

a condition that arises from requiring that the second 
argument of the Hurwitz zeta function satisfies z > 0, 
even for the lowest value of which is obtained for 
u = 0. Furthermore, for the large T expansion to be 
valid, one also requires 

qB/T'^ < 1. (14) 

The conditions expressed in Eqs. m and (ITT)) provide 
the limits of applicability of the high temperature expan¬ 
sion of the effective potential in Eq. m- 


III. THERMO-MAGNETIC COUPLINGS 

We now compute the one-loop correction to the cou¬ 
pling A, including thermal and magnetic effects. Fig- 
ure[T]shows the Feynman diagrams that contribute to this 
correction. Columns (a), (b), (c), (d), (e) and (f) con¬ 
tribute, respectively, to the correction to the cr^, (tt®)^, 
(7r+)^(7r“)^, CT^7r+7r“, (7r'^)^7r+7r“ and tT^(7r°)^ terms of 
the interaction Lagrangian in Eq. 0, respectively. Since 
each of these corrections lead to the same result, we con¬ 
centrate on the diagrams in column (a). Each of the 
three diagrams involves two propagators of the same bo¬ 
son. For the first two diagrams the intermediate bosons 
are neutral and for the third one the intermediate bosons 
are charged. 


The explicit computation has been performed in Ref. . the charged boson Schwinger proper-time propagator for 
This involves the use of the weak field expansion of 
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FIG. 1: One-loop Feynman diagrams that contribute to the thermal and magnetic correction to the coupling A. The dashed 
line denotes the charged pion, the continuous line is the sigma and the double line represents the neutral pion. 


the intermediate charged particles and its corresponding 
qB —>■ 0 limit for the intermediate neutral bosons. Since 
the analysis of Ref. was carried out in the very high 
temperature case, the four-boson vertex correction was 
evaluated in the infrared limit, namely Pi = (0,p —>■ 0), 
where Pi are the momenta of each of the four exter¬ 
nal legs. However, for the present analysis, where T is 
close to Tc a more appropriate treatment is to evaluate 
the vertex function at the typical energy of the external 
particles. This corresponds to the static limit, namely 
Pi = (n,p = 0), where H is given by Eq. (TT^ and rep¬ 
resents the purely thermal (and density) component of 
the boson mass. In the appendix, we explicitly repro¬ 
duce such calculation, evaluating the vertex function in 
the static limit. Notice that Agg- depends on v through 
the dependence on the boson masses. We further con¬ 
sider the approximation where Agg is evaluated at u = 0 
since we are pursuing the effect on the critical temper¬ 
ature, which is the temperature where the curvature of 
the effective potential at = 0 vanishes. 

Next we turn to the calculation of the thermo-magnetic 
correction of the coupling g. Figure [5] shows the Feyn¬ 
man diagrams that contribute to this correction. We are 
interested in computing an effective value for this cou¬ 
pling, ^gg, as well as for n = 0, in the same manner done 
for Agg. Columns (a), (b) and (c) contribute to the cor¬ 
rection to the quark-cr, quark-7r° and quark-7r^ terms of 
the interaction Lagrangian of Eq. respectively. Since 
each of these corrections leads to the same result, we con¬ 
centrate on the diagrams in column (a). Notice that for 
w = 0 in Eq. ®, the masses of the tr and the 7r° become 
degenerate. Hence, the middle and bottom diagrams in 
column (a) of Fig.[2]cancel out, since they contribute with 
opposite signs. This also happens with the two bottom 
diagrams in columns (b) and (c). 

The explicit computation is carried out in the weak 
field limit of the charged boson and fermion Schwinger 
proper-time propagators. The calculation was done in 
Ref. [13 considering only terms of O (qB)'^ after taking 
the average over spins. We have repeated the calcula¬ 
tion and realized that the contribution O (qB) does not 


vanish. This is because such spin average is not needed, 
as we are not considering a decay process, but rather a 
vertex function. We have also evaluated in the static 
limit Pi = = 0) where Pi are the momenta of 

the quark and antiquark and mj'®” is the fermion ther¬ 
mal (and density)-dependent mass, which we compute in 
the appendix and whose square is given explicitly by 


/ them\2 I 

{mf ) = g T [ - - 


1 Lj2(-e^/^) 

3 ~ 


TT^ / 

(15) 


To fix the bare values of the couplings A, g and a ap¬ 
propriate for the description of the phase transition, we 
notice that the boson masses are modified when consid¬ 
ering the thermal effects, since they acquire a thermal 



FIG. 2: One-loop Feynman diagrams that contribute to the 
thermal and magnetic correction to the coupling g. The 
dashed line denotes the charged pion, the continuous line is 
the sigma, the double line represents the neutral pion and the 
continuous line with arrows represents the quarks. 
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FIG. 3: Effective boson coupling Agg- evaluated a the temper¬ 
ature T — 180 MeV with A = 0.4, g = 0.63 as a function of 
the magnetic field strength for different values of g,. 


FIG. 4: Effective boson-fermion coupling ggg evaluated a the 
temperature T = 180 MeV with A = 0.4, g = 0.63 as a 
function of the magnetic field strength for different values of 
g. 


component. For g = 0 they become 


ml{T) 

^liT) 


3Xv^ — + 
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NfN^g^T^ 
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NfN.g^T^ 
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(16) 


At the phase transition, the curvature of the effective 
potential vanishes for u = 0. Since the boson masses are 
proportional to this curvature, these also vanish at = 0. 
In this case, and from any of Eqs. (HSl), we then obtain 
a relation between and the model parameters at the 
critical temperature with g = 0 


a _ jX NfN^ 

TO “ V 2 6 


(17) 


Furthermore, we can fix the value of a by noting from 
Eqs. (HI) that the vacuum boson masses satisfy 


a = 




— 3mJ 
2 


(18) 


TO = 172 MeV [H. Thus, a/T° = 0.77. Equation (d?]) 
provides a relation between A and g. A possible solu¬ 
tion consistent with the above requirements is given by 
A = 0.4, g = 0.63. 

Figure |3] shows the behavior of the effective boson cou¬ 
pling Agg evaluated using T = 180 MeV, as a function 
of magnetic held strength for different values of g. The 
considered temperature is slightly larger that T^. Note 
that Agg is a monotonically decreasing function of qB 
and that the decrease is larger for larger values of g. 

Figure m shows the behavior of ^gg as a function of qB 
evaluated also using T = 180 MeV with A = 0.4, g = 0.63 
for three different values of g. Note that g^g is also a 
monotonically decreasing function of qB. However, the 
decrease is less pronounced for larger values of g. Note 
that the /i-dependence of the effective coupling comes 
from its dependence on 

IV. INVERSE MAGNETIC CATALYSIS AND 
THE EFFECTIVE PHASE DIAGRAM 


Since the effective potential is written as an expan¬ 
sion in powers of a/T we need that this ratio satishes 
a/T < 1. From Eqs. (IT71) and (IT^ the coupling constants 
are proportional to Wo- which, from the above conditions, 
restricts the analysis to considering not too large values 
of nicr as well as not too small values of T°. Since the 
purpose of this work is not to pursue a precise determina¬ 
tion of the couplings but instead to call attention to the 
fact that the proper treatment of screening effects allows 
the linear sigma model to provide solutions for the CEP 
even at finite values of g, we consider small values for 
rria- Given that a is anyhow a broad resonance, in order 
to satisfy the above requirements let us take for defini¬ 
tiveness rria = 300 MeV. namely, close to the two-pion 
threshold. For T° with two light quark flavors we take 


We now study the effect of the thermo-magnetic cor¬ 
rections to the couplings on the critical temperature. For 
a given value of the magnetic field, and for a second order 
phase transition, the critical temperature is determined 
after setting to zero the second derivative of the effective 
potential at r; = 0 in Eq. m- When the phase transition 
becomes first order, the critical temperature is computed 
by determining the temperature where a secondary min¬ 
imum for r; ^ 0 is degenerate with a minimum at r; = 0. 
Figure [S] shows the critical temperature as a function of 
held strength, for different values of g, and for the bare 
values A = 0.4, g = 0.63. Note that in all cases the 
critical temperature is a decreasing function of the held 
strength. 

On the contrary when the calculation is performed 
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FIG. 5: Critical temperature as a function of the magnetic 
field strength evaluated using effective couplings including 
thermo-magnetic corrections with the bare values of the cou¬ 
plings A = 0.4, g = 0.63 for different values of g. Note that in 
all cases the critical temperature is a monotonically decreas¬ 
ing function of the magnetic field strength. 


FIG. 6: Critical temperature as a function of the magnetic 
field strength evaluated without effective couplings and in¬ 
stead with the bare values of the couplings A = 0.4, g = 0.63 
for different values of g. Note that in all cases the critical tem¬ 
perature is a monotonically increasing function of the mag¬ 
netic field strength. 


without including the thermo-magnetic modification to 
the couplings, the critical temperature turns out to be an 
increasing function of the field strength. This is shown in 
Fig. |6] where we plot the critical temperature as a func¬ 
tion of qB for different values of g with the bare values 
of the couplings A = 0.4, g = 0.63. This behavior shows 
that the thermo-magnetic corrections to the couplings 
are crucial to obtain inverse magnetic catalysis. 

We now turn to describe the phase diagram in the 
temperature quark chemical potential plane. Figure [7] 
shows the phase diagram obtained for the bare couplings 
A = 0.4, g = 0.63 for different values of the magnetic 
field strength. The thermo-magnetic corrections enter 
the analysis both in the effective potential of Eq. (HU, 
and in the effective couplings. Notice that as the field in¬ 
tensity increases, the CEP moves toward lower values of 
the critical quark chemical potential, and to larger values 
of the critical temperature and in this case, it reaches the 
T-axis. However, since our analysis is carried out in the 
weak field limit qB/T^ < 1, we can only say that there 
is a tendency for the CEP to eventually reach the T-axis 
large values of the field strength. 

To see the effect of a change of parameters we now ex¬ 
plore the case where the ratio a/T° appearing in Eq. (ED 
changes. We take a/T° = 0.66 which is obtained main¬ 
taining rria = 300 MeV and increasing the value of T° 
to T° = 200 MeV. With this ratio, a possible solution to 
Eq. (1171) for the bare values of the couplings is given by 
A = 0.36, g = 0.51. Figure [5] shows the phase diagram 
thus obtained. Notice that the CEP for qB = 0 happens 
for values of and g'^^^ slightly smaller and larger, 
respectively, than for the corresponding values in Fig. [T] 
When the magnetic held intensity increases the CEP also 
moves toward lower values of the critical quark chemical 
potential and larger values of the critical temperature 


but this time, for the largest value of qB considered, the 
CEP does not reach the T-axis. Nevertheless we observe 
a tendency for the CEP to eventually reach the T-axis 
for larger values of qB that can not be studied within the 
present small held approach. 


V. SUMMARY AND CONCLUSIONS 

In conclusion, we have shown that working in the 
LSMq in the presence of magnetic helds, it is possible 
to obtain values for the couplings that allow to locate 
a CEP that for qB = 0, lays in the region found by 
mathematical extensions of lattice analyses. The analy¬ 
sis is done from the effective potential computed in the 
presence of a weak magnetic held and accounting for the 
plasma screening effects. Since the LSMq does not have 
conhnement, we attribute the CEP’s location to the ad¬ 
equate description of the plasma screening properties. 
Screening is included into the calculation in two manners: 
First, in the effective potential through the boson’s self¬ 
energy and second in the thermo-magnetic corrections of 
the couplings. We have shown that this last correction 
is crucial to obtain inverse magnetic catalysis. To dehne 
the allowed range for the bare coupling constants we ob¬ 
serve that the thermal boson masses vanish at the phase 
transition for g = 0. This condition determines a rela¬ 
tion between the model parameters which can be put in 
quantitative terms from knowledge of and a. The first 
can be obtained from lattice results and the second from 
the vacuum boson masses. Since the model is computed 
in the high temperature limit, we are limited to consider 
ratios of the parameter a/T° a bit off their usual values. 
Nevertheless, the model shows in quantative terms that 
the CEP moves toward lower values of the critical quark 
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FIG. 7: Phase diagram in the temperature quark chemical 
potential plane computed with thermo-magnetic corrections 
to the couplings using the bare values A = 0.4, g = 0.63 corre¬ 
sponding to a/Tc = 0.77 for different values of the magnetic 
field strength. The phase transitions to the left (right) of the 
CEP in each case are of second (first) order. 


FIG. 8: Phase diagram in the temperature quark chemical 
potential plane computed with thermo-magnetic corrections 
to the couplings using the bare values A = 0.36, g = 0.51 cor¬ 
responding to a/Tc = 0.66 for different values of the magnetic 
field strength. The phase transitions to the left (right) of the 
CEP in each case are of second (first) order. 


chemical potential and larger values of the critical tem¬ 
perature as the held intensity increases and that there is 
a tendency for the CEP to eventually reach the T-axis 
for a larger value of the held strength. 

The overall features of the phase diagram can be un¬ 
derstood in general terms when we recall that the mag¬ 
netic held produces a dimension reduction whereby the 
virtual charged particles that make up the vacuum are 
ehectively constrained to occupy Landau levels which, in 
semiclassical terms, implies that their motion is restricted 
to planes. This produces that these particles lay on aver¬ 
age closer to each other. Since as a function of the held 
intensity we have shown that the strength of the inter¬ 
action is reduced, and that this happens no matter how 
weak the external held may be, we infer that a similar 
ehect is taking place in QCD where due to asymptotic 
freedom, the strength of the interaction gets reduced as 
the virtual particles get closer to each other. This weak¬ 
ening of the interaction with proximity between the vir¬ 
tual particles that make up the vacuum should manifest 
itself as well in the weakening of the quark condensate 
with the held strength, as is also observed in lattice QCD 
around the critical temperature. We believe this descrip¬ 
tion will play an important role in the interpretation of 
the lattice QCD results for the behavior of the critical 
temperature and the quark condensate with the held in¬ 
tensity as well as in determining the location of the CEP 
in QCD with and without magnetic helds. 
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Appendix 

A. Thermo-magnetic corrections to the boson 
coupling 


The thermo-magnetic correction to A involves the dia¬ 
grams shown in Fig. [l] column (a). It is only necessary 
to consider the case where the loop is made of charged 
pions (the bottom diagram in column (a) of Fig. [1]), since 
the other contributions can be obtained from this one af¬ 
ter letting H —>■ 0. The calculation is carried out in the 
static limit, i.e. where Pi = (H, p = 0). The explicit 
expression is given by 

/ dP k 
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In the Hard Thermal Loop Approximation (HTL) P 3 y 
P± are small quantities with respect to T, and the same 
occurs with the mass. In this way we find 
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2_Li,2_Li,2 tanh(qBa) ^ 


cosh{qBs) 


Jo cosh(gPT) 


Carrying out the integrals, we obtain 


Jn—oi^ — n) — 


T 


1 


UA*’ I 2 ’ 2 


3 1 , n(r,^) 


IdTT ( 29 P)i /2 
For the non-zero modes (a;„ 7 ^ 0) we find 

dJk 


2qB 


( 21 ) 


.( 22 ) 




n /0 '■ 


1 


(27r)3 

(eBy 


\uj^ + k^ + (uj^ + k^ + 


2(eBfk 


2 j ,2 


(w^ + k^ + m^)'^ J V ~ ‘^n)^ + (p — k)^ + 

{eBf 


((w - uJnY + {p- ky + rrpy 


+ - 


2{eByP 


((w - ujnY + (p - ky + m^yj 
Still in the HTL approximation we find 

n^O'' 

1 {eBy 


(23) 


(2^)3 


uj^ + 2 {lu^ + k'^y 

1 


{eBf 


{uj-LUny + k'^ {{lU - UJny + ky^ J _ 


(24) 


Since we only consider terms up to order O {qBy, we 
have 


Jn^o{Pi') TTii) — T 'y ] 

n# 0 ‘ 

1 


d^k 

W 

{eBy 


{uj^ + k^ + uj'^y {u}^ + k'^ + ^ 2)4 


(25) 


In order to calculate the above integrals we make use of 
dimensional regularization, and of the Mellin summation 
technique [s^, to find 


Joining both contributions we find 
T 1 

J(w = n) = 





In the case of the diagrams involving neutral bosons we 
have 


/(P,;m 2 ) = f -^D{P,-K)D{K) 

= In=o{Pi]mJ) + In^yPi^mj). (28) 

In order to calculate I{Pi;mf), we take limit (qB) —>• 
0 in Eq. ( 122 ]). The limit of the Hurwitz Zeta function 
is not trivial and we we use the following assymptotic 
expansion [s^ 


Cis,y) = " 


y^ ^ ^ B2k T{2k + S-1) 

s - 1 ( 2 fc)! r(s)p 2 fc+«-i ’ '' 


where B 2 k are Bernoulli numbers. This expansion is valid 
for large values of y, which is equivalent to having a small 
value for qB. In our case s = 3/2 and y = ^ + and 
we find 

^/3 1 H \ 2(2gB)i/2 1 {2qBf/'^ 

^[ 2 ’2^^) ^ 7n 16 n5/2 ■ 

(30) 

Using the above expansion in Eq. (EZj) we obtain 


I{u} = U) = - -^ 


T 1 

Stt Vn 

1 r. /( 47 rT )2 


IGtt^ . 


In 


C(3) 


'Vn' 

2'kT 


2 a 2 

2 


-I- 1 — 2^e 


(31) 


B. Thermo-magnetic corrections to the 
fermion-boson coupling 

The determination of the thermo-magnetic correction 
to the coupling g involves the diagram shown in Eig.|2](a). 
We call the one-loop effective vertex T. The calculation 
is done up to order O {qB) 


Jn^o{^ — H) — 


+ C(3) 



1 


In 


(^) 


- 1-1 — 2'^e 


{^py 

10247r6T4 


C(5). 


(26) 


r = -3 

■5-^ f dJk 

+ j^1oS{Pi-K)S{P2-K)^oD{K) 

= -g{l + ST). (32) 
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Concentrating on 

<5r = I 


(Pk 


(27r)3 

X S{P2-K)-i^D{K) 

n f d?k 


755(Pi - K) 


-{fi-m 


(Pi - KY + m? 


*9^7172(7 ■ (-P2 - i^)||)1 r - if 2 - fC ) 

[(P2-P:)2+m2]2 J |_(P2 - P:)2 +to' 

iqBjij2i'J ■ (P 2 - ^)||) 


[(P2-P:)2+m2]2 
In the HTL approximation, we get 

(Pk 




m. 


2 ■ 


(33) 


5r = -9^tY, 


(^)^A(Pi - iF)A(P 2 - K) 


(2^)3 

X A(iF)-iqP^7i72(7-/0||A(Pi-iF) 

X A2(P2-P)A(iF)-^gP7l72(7•i^)||-^ 

X K{Pi-K)K^{P2-K)l^{K) 

= + 6Ttb, 


(34) 


where 


(5r 


TV 


= -sTY, 


(Pk 




(2^)3' 

X A{Pi-K)A{P2-K)A{K), 
is the vacuum + thermal contribution, and where 

(Yk 


(35) 




TB 


= -sTY 


X A{Pi-K)AY{P2-K)A{K) 
- iqB"fi"f2{l ■ K)\\^A{Pi - K) 

X A^{P2- K)A{K) 


-iqBfC^i'j2{'l ■ K)\\ 


(36) 


is the thermo-magnetic contribution. We now consider 
i 5 r 7 ’y in the HTL approximation 

^Ftv = J ^^^A(Pi-iF) 

X A{P2-K)A{K) 

= 8^ 

Next, we concentrate on the last two terms of Eq. (El, 
which make up the thermo-magnetic contribution STtb- 
First, we recall 


(38) 


anti-commutes with the other gamma matrices. We in¬ 
troduce the decomposition 

7172-^11 = 75 [(A'- 6 )^"-(i^- m)^ , (39) 

where we have introduced the four-vectors 


= ( 1 , 0 , 0 , 0 ) 

= (0,0,0,1). (40) 

We stress that in the HTL approximation, Pi y P 2 are 
small quantities that can be considered of the same or¬ 
der. In this way the thermo-magnetic contribution can 
be written as 

STr^ = -g‘TY f - if) 

X A 2 (P 2 - K)A{K)[-2iqB-f5{T/iK ■ b)) 

- m-uW- (41) 

We define 

r 93 u 

G(Pi,P2) = tJ2J J^MPi-K)AYP2-K) 

X A{K)MK-b))-m-uW: (42) 

to obtain 


STtb = 2ig\qB)j5GiPi, P 2 ), (43) 

where G(Pi, P 2 ) can be expressed in terms of the tensor 
Jai (a = 1 ,... 4, i = 3,4) given by 



X A^{K)A{Pi-K)A{P2-K). (44) 


In order to calculate the sum over Matsubara frequencies 
we use 


lo 


El 


E 2 


a2(p:)A(Pi - p:)A(P2 - a:) 

n 

Tj2m)A{P, - K)A{P2 - K) 

^ ' n 



T^w„A2(/F)A(Pi - K)A{P2 - K) 

n 

tY,^Mk)A{P, - K)A{P2 - K) 
T^ujlA^{K)A{Pi - K)A{P2 - K) 

n 

(-^) EE^nA(P)A(Pi - K)A{P2 - K) 



75 = 74717273 


(45) 
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where Xq, Xi, X 2 are given by 


E SS 1 S 2 1 

SEE1E2 — ^2) — ^ 1-^1 H“ S2E2 

S,Si,S2 

' l-fisE) + fisiE,) l-J{sE) + fis2E2) 
iuJi — sE — siEi iuj 2 — sE — S 2 E 2 

E S 1 S 2 E 1 

SEE1E2 i(^i — ^2) — ^lEi H- S2E2 

fi «T fin ^ ' 


~ 1 - fi^sE) + /(si^i) 1 - f{sE) + f{s 2 E 2 ) 

iuJi — sE — siEi iuj 2 — sE — S 2 E 2 

E SSiS2E‘^ 1 

SEE1E2 — <^ 2 ) — siEi + S2E2 

,Sl,S2 

' 1 - fjsE) + fjsiEi) 1 - f{sE) + f{s 2 E 2 ) ' 
iuii — sE — siEi iuj 2 — sE — S 2 E 2 

(46) 


The leading temperature behavior is obtained from the 
terms with s ^ —si = —S 2 - We consider in detail the 
calculation of Xq for those terms and make the approx¬ 
imation where f{Ei) ~ f{E 2 ) — f{E), namely, that the 
Bose-Einstein distribution depends on E = -b 
and thus on the quark mass. This approximation allows 
to find the leading temperature behavior for m —> 0, 
since it amounts to keep the quark mass as an infrared 
regulator. Also, using that E^ ~ k — pi ■ i = 1,2, we 
find 


1 [f{E) + f{E) 

~W E 

1 

{iuji +pi ■ k){ioj2 + P2 ■ k) 

1 

^ 

(iwi — Pi ■ k){iuj2 — P2- k) 




where we have set Ei = E 2 = k in the denominator of 
the first fraction. Similarly 


Ai ~ - 


f{E) + f{E) 


8k 


E 


1 


[ {iuJi +pi ■ k){ioj2 +P2-k) 

_d_.1 

{liOl - Pi ■ k){lUJ2 -p2-k) I 




fi.E) + f{E) 


E 


1 


{ioji +pi ■ k){ioj2 +P2- k) 

1 1 


(48) 


(iuJi — Pi ■ k)(iuj2 — P2 ■ k) J 

Using Eqs. (ITTl) and (H5)) in Eqs. (IT5|) and (HTl) . we find 

dx 


ai — 


—f -—1 r 

87r2 ^ dy^J Jq 




f(y/xE+^) -f f{y/FFE^) 
f dn kaki 


4^ {Pi-k){P 2 -ky 


(49) 


where we have defined x = k/T, y = m/T, K = {—i, k), 
Pi = {—uji^pi) and P 2 = {—ui 2 ,P 2 )- The integrals over x 
can be expressed in terms of the well known functions [s^ 


hn{y) = 
fn{y) = 


1 


1 


r(n) 

1 


3 yjx^ Py"^ gy/x'^+y'^ _ ^ 

dxx^-^ 1 


r(n) Jo yjx^ - 1 - 2/2 g^a: 2 -|-y 2 ’ 

which satisfy the differential equations 

^^n+l ^n—1 


(50) 


dy"^ 

dfn+l 

dy^ 


2n 

fn-l 

2n ’ 


(51) 


therefore 


3OLl - 


1 


[hliy) + fi{y)] 
4 ^ {Pi ■ k){P2 ■ k) ■ 


IGtt^ 
' dVL 


(52) 


Using the high temperature expansions for hi{y) and 

h{y) 0 

= i; + )‘"(s)+r* + '" 

h(y) = -)l”(f)-)7* + ..., 


(53) 
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and keeping the leading terms, we obtain 


^Oii 


Hence 


IGtt^ 


ln(2) - 

, 2 yn. 


{P, . K){P^ . K)^ 


STtb = — 




IGtt^ 

dn [{^{k ■ b)) - ^{K ■ u)]k 
4^ (Pi • k){P2 ■ k) 


ln(2) - 

. 2 Vn. 


(55) 


In order to consider the thermo-magnetic dependence 
of the fermion-boson coupling we consider explicitly the 
quantities appering on the r.h.s of Eq. (1551) 


d^Qj(Pi, P 2 ) = 




K„K,. 


47r (P, . K)(P2 ■ K) 


(56) 


For simplicity we choose a configuration where the mo¬ 
menta Pi and p 2 form a relative angle 0i2 = tt. This con¬ 
figuration corresponds, for instance, to a thermal gluon 
decaying into a quark-antiquark pair in the center of mass 
system, and is therefore general enough. Consider first 
JikPl,P2) 


J44(Pi,P2) = — 


1 


1 

2 iuJiP2 + ioJ2Pi 

f dx\ - - - 

_i ( iuji + pix 

1 1 

2 iuJiP2 + iw2Pi 

iuji +pi 


P2 


In 


ILOi - Pi 


In 


1002 — P 2 X 


/ iU!2 + P2 
\iUJ2 -P2 


(57) 


We now perform the analytic continuation to Minkowski 
space iooi ^2 —t Poi ,02 [k —>■ (—l,fc)], and consider the 
scenario where poi = P 02 = Po and pi = p 2 = p, leading 
to 


^44 -t ^00 = 


2poP 


Po-P 


(58) 


Furthermore, we consider the static limit where the 
quarks are almost at rest, namely p —>■ 0, to find 


p-iO 1 

Joo —' — 2 - 

pI 


(59) 


Now we consider JssiPi, P 2 ) in the same momenta con¬ 
figuration 

J33{Pi,P2) = 


2 iojip2 + ioJ2Pi 
X / dx x'^ 


1-1 


Pi 


P2 


ioui -I- piX ioj2 — P 2 X 


1 


i0OiP2 + i0O2Pl 
iuji 

2pi 


Pi L 
i0J2 
P2 


lOJi ( iu!i+ Pi 

1 — 1 :— In 

zwi - Pi 


[l-^ln 

f iW2 + P2 ^ 

2p2 

V 1002 - P2 / 


(60) 


After analytical continuation to Minkowski space and in 
the same scenario where poi = po 2 = Po and pi = p 2 = p, 
we obtain 


J 33 = — - 


1 

p2 



( Po+p\ 

L 2p 

O 

1 


In the limit where p —>■ 0 this gives 

p-lO 1 

J 33 —t 

m 

In this same limit, p —>■ 0, we find 


STtb = 2p^(J33 -|- J44)q'^ ■ B 


(61) 


(62) 




167r^ 


= 2-7^ ( :t-2 


4 

3Po 




167r^ 


■gS • P, (63) 


where 


T, ■ B = 17172 P. 


(64) 


By taking both contributions into account, STtv + dTxBi 
we can find F. Considering the contributions of the up- 
and down-quarks explicitly, and taking Pq —?► mj we fi¬ 
nally obtain the correction to the coupling g. 


Self = 5(1 + 5^ (5TB + 5Tv))> 


(65) 


where 


ttT \ 


{qu + qd)B f f, 

Stt^ [smp) i ^ ^ 2Vn; 

9TV = (in (^) + f - ) - ln(2»)) ■ (66) 
C. Fermion thermal and density dependent mass 


This calculation involves the three diagrams shown in 
Fig.O We only consider the first one, since the computa¬ 
tion of the other two diagrams is completely equivalent. 
We call this diagram 
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For the angular part we notice that the integral is sym¬ 
metric under the transformation 


implying 


k —^ — k‘^ dfl —y dfl, 


(73) 


dfl 1 

47r — fc .' 


dfl 1 

47r iui + k ■ p 


(74) 


FIG. 9: Diagrams contributing to the calculation of the 
fermion thermal mass. 


=-9^ I - K)MK). 

In the HTL approximation this gives 


where 


A(g) 

A(Q) 


1 

LOn + k^ 

1 

Uln + k'^ 


(67) 


( 68 ) 


(69) 


The determination of involves two kinds of Matsubara 
sums 


= tJ2Mk)Mp-k) 

n 

= Tj2^nA{K)AiP-K). (70) 


The first case refers to i = 1,2,3 and the second case to 
1 = 4. The first kind of contribution is 


= ^E 


d^K 


1 r dn 

27r^ J Att 

f{k) + f{k + p) 


^K,A{K)A{P - K) = 


Ki / dkK 


f{k) + f{k - p) 


iio + k ■ p 


iuj — k ■ p 


(71) 


Therefore, we finally find 


E(Ti = i^2,3 


Stt^ 


(^)-T2E^2(-e'^/^) 




dn 1 

47r ioj + k ■ p 


• (75) 


The second contribution is 


= ^E 


d^K 


dkK 


(27r) 

i f dQ 

Stt^ J drr 

f{k) + f{k + p) 
ioj — k ■ p 


^uJnA{K)A{P - K) = 


f{k) + f{k - p) 


iu! + k ■ p 


(76) 


Carrying out the radial and the angular integrals, in a 
completely analogous way as before, we have 


= ^E 


d^K 

^2rp2 


ujnA{K)A{P - K) 




- r2Lf2(-e”^/^) 


dn 


47r ioj + k ■ p 


(77) 


The sum of the contribution of the different diagrams fi¬ 
nally yields the thermal and density dependent correction 
to the fermion mass 


We have to deal with both the radial as well as the an¬ 
gular part. For the radial contribution we find 




dK K[f{k) + f{k ± p)] = 




1 Li2(-e^/^) 

3 TT^ 


Li2{-e ^/^) 


(72) 


(78) 
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